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^ . Abstract 
> . 

^ ' Invariant theory as a study of properties of a polynomials under translational trans- 

formations is developed. Class of polynomials with congruent set of eigenvalues is in- 
troduced. Evolution equations for eigenvalues and coefficients remaining the polynomial 
^ ^ within proper class of polynomials are formulated. The connection with equations for 

■ hyper-elliptic Weierstrass and hyper-elliptic Jacobian functions is found. Algorithm of 

. calculation of eigenvalues of the polynomials based on these evolution equations is elab- 

orated. Elements of the generalized dynamics with n-order characteristic polynomials is 
built. 

in 
o 

00 ! Keywords: polynomials, elliptic functions, algorithm, evolution equation, dynamics, rela- 

^ ' tivistic mechanics. 
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j_l 1 Introduction 



The problem of expressing of eigenvalues of the polynomials as a certain functions of the coeffi- 
cients is one of the oldest mathematical problems. The question on possibility, or impossibility, 
to express the eigenvalues of the polynomials through coefficients on making use of radicals had 
been exhaustively answered by E.Galois and H.Abel that the polynomial higher than fourth 
order, in general, does not admit a presentation of solutions via radicals [3]. In spite of this 
rigorous mathematical theory mathematicians remained to believe that the eigenvalues of the 
polynomials can be expressed in analytical way as certain functions of the coefficients [2], [7]. 
Hermite was first who found a very elegant expression of eigenvalues of the quintic equation by 
modular functions [6] . The theory of elliptic functions originally was related with the problem 
of finding of eigenvalues of the cubic polynomial. In fact, the Weierstrass elliptic functions at 
the periods are equal to eigenvalues of the cubic equation [12]. It is clear, however, that for a 
search of analytical solutions of the n > 5-degree polynomials one needs of mathematical tools 
beyond the elliptic functions. In this context as a hopeful tools one may take the theories of 
hyper-elliptic functions [9] and multi-complex algebras [8], [13]. 
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The main purpose of the present paper is to construct evolution equations for eigenvalues 
and coefficients of polynomials. We search such kind of evolution equations which remain the 
original polynomial within certain class of polynomials. Therefore first of all we define the 
set of invariants and classify the polynomials with respect to the obtained set of invariants. 
The polynomials from the same class possess with congruent set of eigenvalues, i.e. inside a 
given class the eigenvalues of one polynomial are obtained by simultaneous translations of the 
eigenvalues of the other polynomial. The algorithm of calculation of the eigenvalues is based 
on the evolution process reducing the number of the coefficients of initial polynomial. During 
of the evolution the initial polynomial will transformed into the other polynomial remaining 
within the frames of a given class. The evolution is directed in such a way that the final 
polynomial will possess with one trivial solution. The coefficients of the finial polynomial are 
solutions of the Cauchy problem for ordinary differential equations where the coefficients of the 
polynomial serve as initial data. As soon as the Cauchy problem is resolved, the eigenvalues 
of initial polynomial are found simply by certain set of translations of the eigenvalues of the 
final polynomial. It is shown, if the coefficients obey to equations for Weierstrass hyper-elliptic 
functions then the eigenvalues will obey equations for hype-elliptic Jacobi functions. 

The present method has been procreated in the process of construction of generalized elec- 
trodynamics of n-th order (see, Refs.[14],[15],[16]). In this theory an important role plays the 
fact that the mapping between inner- and outer-momenta is built as a mapping between coef- 
ficients and eigenvalues of the characteristic n-degree polynomial. The generalized dynamics is 
based on dynamic equations of motion previously developed for the polynomials. As a starting 
platform of the construction served the elements of the relativistic dynamics closely related 
with quadratic polynomial. 

Besides the Introduction the paper contains the following sections. 

In Section 2, the equations of evolution for the coefficients of n-degree polynomial are for- 
mulated. In Section 3, the Algorithm of finding of eigenvahics of the corresponding polynomial 
is built. In Section 4, some peculiarities of the cubic equation is explored. In Section 5, the 
elements of the relativistic dynamics based on quadratic characteristic polynomial is presented. 
In Section 6 we give an account of a sketch of the generalized dynamics related with n-degree 
characteristic polynomial. 

2 Evolution equations for eigenvalues and coefficients of 
n-degree polynomial 

If F is a field and qi, ...,qn are algebraically independent over F, the polynomial 

p(X) = fl{X-q,) 

1=1 

is referred to as generic polynomial over F of degree n. The polynomial equation of n-degree 
over field F is written in the form 

p{X) := X" + J2H'in -k+ l)PfcX"-'= + (-)"P2 = 0, (2.1) 

k=l 

where the coefficients Pk G F{qi, ...,qn)- In this paper we shall restrict our attention only to 
polynomials with real coefficients and with simple roots. Moreover the last coefficient is 
essentially positive. The signs at the coefficients in Eq.(2.1) are changed from term to term 
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which allows in Vieta formulae to keep only the positive signs. The expressions at the coefficients 
are included purely for convenience, and have no real bearing on the theory. The mapping from 
the set of eigenvalues onto the set of coefficients is given by Vieta formulae 

n n k 

(a) nPi = ^g,, ih) P^ = ^qu (c) = E ^Qr,, (2.2) 

1=1 1=1 \<ri<...<rk<n j=\ 

here is called the elementary symmetric k-degree polynomials of the eigenvalues. The number 
of monomials inside k-th elementary polynomial is equal to binomial coefficient: 

^'={t) = k\{n-k)V ^^-^^ 

Since the roots of the generic polynomial p{X) are algebraically independent, this polynomial 
is, in some sense, the most general polynomial possible. 

In p{X) replace X with X = Y + Pi. This transformation will eliminate the (n — l)-degree 
term, resulting in a polynomial of the form 

n-l 

r(r) := y" + Y.H'RkY''-'' + i-Y'Ro = 0. (2.4) 

k=2 

The polynomials p{X) and r{Y) have the same splitting field and hence the same Galois group. 
Let E be the splitting field of r{Y), let y^. A; = 1, .., n be its roots in E and G — Gf{E) be its 
Galois group. 
Lemma 2.1 

The coefficients Rk, Rq of Eq.(2.4) are invariants with respect to simultaneous translations 
of the eigenvalues of Eq. (2.1) 

The Proof of the statement it follows from formula Y — X — Pi which in terms of the 
eigenvalues is expressed as follows 

Vk^qk — ^qi^ -Y^{qk- qi)- (2.5) 

1=1 i^k 

It is seen that the eigenvalues of Eq.(2.4) are represented by differences between the eigenvalues 
of Eq.(2.1), hence, they arc invariants with respect to the simultaneous translations. Since the 
coefficients RQ,Rk,k = 2, ...,n — 1, are sum of uniform monomials of y^. A; = l,...,n they 
have same feature, namely, they are invariants with respect to simultaneous translations of the 
eigenvalues of Eq.(2.1), too. 
End of Proof. 

The polynomial r{Y) we denominate as invariant polynomial (with respect to translations 
of the roots of (2.1)). 

The main task of this section is to introduce evolution equations for the coefficients of 
Eq.(2.1) which remain invariant the coefficients of Eq.(2.4) Rk. This result is given by the 
following 

Theorem 2.2 

Let qk,k = l,...,n be set of eigenvalues of polynomial equation of n- degree (2.1). Let the 
differentials of all eigenvalues are equal to each other 

dqi = dq2 ^ . . . ^ dqk = . . . = dqn- (2.6) 
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Then differentials of the coefficients satisfy the following system of equations: 



2Pn-idPi = dP^, 



(2.7) 



dPn-k — 



{k + 2)P^_k-idPi, k^l,...,n-2; 



(2.8) 



Proof. 



Notice that from (2.2a) it follows that 



dqk = dPi, k = 1, ...,n. 



Coefficients of the polynomial (2.1) are symmetric forms of the eigenvalues where k-th coefficient 
{n — k + l)Pk consists of monomials of /c-degree. Thus, the derivative of this coefficient 
{n — k + l)dPk contains kC^ monomials and, since the derivatives all eigenvalues are equal to 
each other, so, the derivative of {n — k + l)Pk is proportional to XdPi where the coefficient of 
proportionality consists of kC^ {k — l)-degree symmetric monomials. On the other hand, the 
symmetric polynomial of {k — l)-th degree can expressed only by C^~^ symmetric monomials of 
(/c — l)-degree. This means the expression for A consists of kC^JC^~^ = n—k+1 same symmetric 
polynomials of {k — l)-degree which are equal to the {k — l)-th coefficient {n — k + 2)Pk-i. The 
result is expressed as follows 



which completes the system of differential equations for the coefficients of Eq.(2.1). 
End of Proof. 

The following Lemma demonstrate an important role of the invariant polynomial in the 
evolution process. 
Lemma 2.3 

The first integrals of evolution equations (2.8) are given by coefficients of invariant polyno- 
mial (24). 
Proof 

In fact, the roots of Eq.(2.4), i/^, k = l,...,n, according to formulae (2.5), are the first 
integrals of evolution equations for eigenvalues (2.6). Equations (2.7)-(2.9) are consequences 
of Eqs.(2.6), hence coefficients Rk as algebraic functions of the solutions of evolution equations 
are first integrals of Eqs.(2.7)-(2.9). 

End of Proof 

Inversely, the use formula Y = X — Pi in Eq.(2.4) has to transform invariant equation 
into equation (2.1). Substitute X — Pi instead of Y and gather together powers of X and 
the coefficients of the obtained polynomial compare with coefficients of Eq.(2.1). We shall see 
that the coefficients P^, k = 1, ...,n — 1 now arc expressed as /c-degree polynomials of Pi with 
coefficients consisting of Rk, k = 2, ...,n — 1. Especially notice, the invariant Rq is defined by 
n-degree polynomial of Pi with coefficients built from Rk- The first task is to find an explicit 
form of this polynomial. The general form of those polynomial can be presented as follows 



{n-k + l)dPk = {n - k + l){n - k + 2)Pfc-i dPi, k = 2,3, n-1. 



Differentiation of the last coefficient gives 



dP^ = 2P„_idP] 



(2.9) 



n-1 



P" + ^ P" ^/fc(Pi, Rk) + Ro — P^- 



k=2 
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Now, the task is to find an explicit form of the function f^. With that purpose explore firstly 
the case = 0. In this case one of the solutions of Eq.(2.1) is equal to zero. Then the 
corresponding solution of the invariant polynomial is 

yi = -Pi(p2 = 0) = -Pi(0). 

Hence (-Pi(O)) will satisfy (2.4). By replacing Y with (-Pi(O)) in Eq.(2.4) we come to the 
following equation for Pi. 

n-1 

pr(o) + E RkPr\o) + i?o = 0. (2.10) 

k=2 

Notice, here the signs at all coefficients are positive. Secondly, suppose that we have changed 
coefficients of Eq.(2.1) from the set { P^, 7^ } to the set { P^, P^ = } obeying evolution 
equations (2.8). This way of evolution provides the polynomial with P'^ — with the same 
invariants as the original one. Hence, when P'^ the coefficients of Eq.(2.10) will not change, 
but now this polynomial will be equal to P^: 

n-1 

P^ + Y.RkPi~'' + Ro = P^. (2.11) 

k=2 

Notice, the situation is somewhat similar the conventional Classical Invariant Theory of Poly- 
nomials [10]. If classical invariant theory is a study of properties of a polynomial p{x) that 
are unchanged under fractional linear transformations, within the framework of the present 
approach we study properties of polynomials under translational transformations. 
Theorem 2.4 

Let coefficients of Eq.(2.1) obey evolution equations (2.7)-(2.9). Then 

(^+l)^n-.= ^^P^. (2.12) 

Proof. 

Differentiate equation (2.11) by taking into account that Pq; Rki k = 1, n—1 are constants. 
We get 

n-2 

dPi ( nPf-^ + ^{n- k)RkPr''~^ + Rn-i ) = dP^. (2.13) 

k=2 

Compare this equation with equation (2.9). It is seen, the expression inside brackets at dPi in 
(2.13) is nothing else than the coefficient 2P„_i expressed as a polynomial of Pi with invariant 
coefficients: ^ 

2P„_i = nPr' + - k)RkPr'~' + Rn-1- (2.14) 

k=2 

Hence, 



2Pn-l = ip^. 



dPi 

Next, differentiate (2.14) with respect to Pi, we obtain 

dPn-i = 3P„_2 dPi, (2.15) 
where 3P„_2 in fact is the next coefficient of Eq.(2.1): 

3P„-. = ^^t^PT-' + "£ " fl./T" + Rn-. (2.16) 

^ k=2 ^ 
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Hence, 

At the next step we shall obtain 

dPn-2 = ^Pn-3 dPi, (2.17) 

where the expression at the differential dPi we denoted by 4P„_3 because this expression indeed 
is {n — 3)-th coefficient of Eq.(2.1): 

4P„_, = n(»-0(B-2) ^_3 ^ g („-.)(„-.- l)(n-.-2) ^^^„_..3 ^ 
Hence, 

From these formulae by induction one may easily establish that the general formula for l-th 
coefficient Pn-i is 

(' + 1)^"-' = ( n ) ^' + ( „ - + ^»-' = S Jv""- 

End of Proof. 

This theorem has some interesting consequences. 
Corollciry 2.5 

The following representation for polynomial p{X) holds true 

exp(-X^)p2 ^ 0. (2.19) 

Proof 

The Euler operator, generator of translation, is represented by following expansion 

exp(-X— ) = 1 - X— + X^l^ + ... + X^l^ + ... . (2.20) 

By differentiating Eq.(2.11) n-times we get 

-^P2 = n\. (2.21) 
dPf ^ ' 

Hence the sum in (2.20) is completed by this term. Thus, 

exp(-X^)p2=p(X). (2.22) 

Here the variable X means one of the roots of Eq.(2.1) in quality of which let us take Qn- 
The derivative with respect to Pi due to (2.6) can be expressed by the following sum 

J_-d^d_-d_ 

dPi tidPidqk ti^Qk' ^ ' ' 
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Then, 

d d d 

exp(-5„— )p2 = II exp(-?„— ) exp(-g„— ). (2.24) 

Now, take into account Vieta formula for given by (2.2). The Euler operator (2.24) will 
translate each root by g„. The last operator in (2.24) acts only upon n-th root resulting 
?n - ?n = 0. Hence, 

exp(-X^)p2=p(x) = 0. 

End of Proof. 



3 Algorithm of finding the eigenvalues of n-degree poly- 
nomials 

The main idea of the present algorithm is to reduce the problem of solution of n-degree polyno- 
mial equation into the problem of solution of n — 1 degree polynomial equation. The evolution 
from n-degree polynomial up to (n — l)-degree polynomial is fulfilled in such a way that remains 
invariant coefficients Rk- Hence, initial and final polynomials of this evolution will possess with 
congruent set of eigenvalues, so that the solutions of the former can be obtained from the 
solutions of the latter simply on making use of transformations of translation. 

One of the ways to reduce a degree of the polynomial is achieved by tending P^ to zero. For 
that purpose we must to use P^ as an evolution parameter of the evolution equations (2.8) with 
final goal to find the coefficients of Eq.(2.1) at the point P^ — 0. This evolution will remain 
invariant the coefficients of Eq.(2.4), hence from solutions of the polynomial with = we 
may come to the solutions of the original equation simply by translation of the set of solutions 
with new coefficients. 

Denote x = P^. Re-write Eq.(2.8) with respect to x. We get 

2Pn-idPi^dx, , 
dPn-k = {k + 2)Pn-k-idPu /c= l,...,n-3; (3.1) 

P = nP.dP,. 

ax 

This a well-known Cauchy problem with initial data Pk{x — P^), k — 1,2,3, ...,n — 1. 
The variable x run from x — P"^ till x — 0. These equations usually are resolved by using the 
celebrated Cauchy-Lipschitz method of calculation [4] . This procedure is fulfilled by dividing 
the interval {xq, 0) into parts: 

Axo = - Xq, Axi = Xi+i - Xi, Aa;„_i = x - x^-i, (3.2) 

where Xi < Xj+i, xjy = 0. 

In this way the continuous evolution process is transformed into discrete process consisting 
of N steps. In the last step of we come to n-degree polynomial free of the last coefficient: 

(1) n-l (1) (1) 

X'' + J2{-)''{n-k + l)PkX'' ^0. (3.4) 

k=l 
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One of the solutions is trivial, excluding this solution we come to the polynomial of (n — 1)- 
degree: 

(1) n-l (1) (1) (1) 

+ J2i-)\ri -k + l)PfeX'= + 2P„_i = 0. (3.5) 

k=l 

Let us mention that Eq.(3.4) possesses with same invariants as original one, i.e. Eq.(2.1). If 
one will use the numerical methods of solution then the relationships for invariants are satisfied 
within given accuracy of calculations. Suppose that the solutions of {n — l)-degree equation 

(1) (1) 
(3.5) Q k, k = 1, ...,n — 1 are known. Complete this set of solutions by ? „ = 0. Then the 

solutions of original equation one may find simply by the following translations 

(1) (1) 

Qk^ Qk + Pi- Pi, k^l, n. (3.6) 

If solutions of the (n — l)-degree equation still are unknown, then one may continue to apply 
this algorithm again in order to reduce the problem of solution of (n — l)-degree equation to 
the problem of solution (n — 2)-degree polynomial equation. This process can be continued 
up till quadratic or linear equation. At each step of iteration one will find an information on 

(r) 

coefficient Pi. At the r-th iteration one deals with (n — r)-degree equation in the form 

(r) n-r-1 (r) (r) (r) 

fe=l 

(r) 

At this stage the norm of the coefficient Pn~r will fulfil the role of evolution parameter of the 
next evolution process. As soon as the solutions of the lowest polynomial is found, the inverse 
process of iterations will consists only translations of the known set of eigenvalues with known 
parameter of translation. Let the last step of iteration be linear equation from which we find 
only one solution, 

(„-i) ('^-i) 
gi = n Pi . 

Then the solutions of the original equation (2.1) are found as a result of the following set of 
translations 

(n-l) " \ (n-r) (n-r+1) 

qi^nPi +nY^-{ Pi - Pi ), (3.8a) 

r=2 ^ 

" X (n-r) (n-r+1) 

g. = E-( ^1 - ^1 , s = 2,3,...,n; (3.86) 

r=s ^ 

(0) 

where Pi = Pi . 



4 Relativistic Lorentz-force equations and related quadratic 
polynomial 

Inter-relation between the evolution of the eigenvalues an the coefficients of the polynomials 
with dynamic equations for physical systems we can observe already at the level of quadratic 
polynomial. In aim of this section is to demonstrate as the quadratic polynomial is related 
with relativistic Lorentz-force equations. In the sequel, we shall use this example as a starting 
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platform to pass to case of higher degree polynomials and related generalized dynamics. Let 
us start with generic quadratic polynomial 

-2Po X + ^0, (4.1) 

with real coefficients Pq > P^, and real eigenvalues This polynomial closely related with 

the relativistic dynamics [14], [5]. In order to demonstrate this fact let us start from Lorentz- 
force equations for a charged particle inside external electromagnetic fields E, B written with 
respect to proper time r ^ 

— = (EPo + [PxB]), ^ = (E-P), (4.2) 
dr ^ dt ^ , , , 

Consider projection of equations (4.2) on direction of the motion defined by unit vector n — 
P/P: 

^£.(E.n)P„,f ^(E.n)P, (4.4) 

Then one deals only with the lengths of the momenta Pq, P- Simplify Eqs.(4.4) by introducing 
new evolution parameter: 

where 

^ = (E-n). (4.6) 
The first constant of motion is easily found 

p2 _ p2 ^ M^. (4.7) 

Here is a constant of motion. Conventionally, this constant is interpreted as a square of 
inertial mass [11]. Inside stationary potential field, when E = —'W{r), dynamic equations 
imply the other integral of motion, the energy, 

So = Po + V{r). (4.8) 

At the rest state where P — one obtains Po(P = 0) = M. The relativistic mechanics deals 
with two kinds of the energy, namely, 

pj^Po-M, pI^Po + M. (4.9) 

From formulae (4.7) and (4.9) it follows 

P^PiP2, Po^I{pI+pI), M^^{pI-pI). (4.10) 

Notice, the first two formulae of (4.10) mean Vieta formulae for quadratic polynomial equation 
(4.1). Substitution X = Y + Pq in (4.1) leads to invariant equation 

= P^ - P2 = M^, (4.11) 



^Hereafter for the sake of simplicity we omit all parameters like charge, mass, light-velocity and other 
parameters regularizing physical dimensions taking them equal to unit. 
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the invariant coefficient of which is equal to the invariant of physical motion. Quadratic equation 
for Po is given by 

Differentiating this equation we derive evolution equations for the coefficients, which in the case 
of quadratic polynomial, of course, is a simple task: 

2Po P = ^Po = P (4.13) 

as as 

Solutions of these evolution equations are given by hyperbolic cosine-sine functions 

P ^ M sinh{s), Pq^ M cosh{s). (4.14) 

Then the eigenvalues of Eq.(5.1) are expressed by hyperbohc cosine-sine functions of one-half 
argument 

pl^V2Msinh{^), pj ^ V2M cosh{^) . (4.15) 

Now, consider so-called the effective potential representation. With this purpose let us come 
back to the equations written with respect to proper-time r. Then from the second equation 
of (4.13) it follows 

£o = Po + Vir). 

Further, replace Pq by £o — V{r) in the first of equation of (4.13), this gives 

= _vy(r) {So - V{r) = -VW{r, So), (4.16) 

where the effective potential is defined by 

W{rSo)=SoV{r)-^V\r). (4.17) 

Here, the relativistic Lorentz-force equation written in the form of Newtonian equation with 
the effective potential W. From this equation it follows the Newtonian form of the energy 

S^Ip' + W{r, So) = liS'o - M'). (4.18) 



5 Cubic polynomial equation and related dynamics 

In the previous section we have demonstrated as the relativistic dynamics is related with the 
evolution of quadratic polynomial. This example provides us with appropriate tool in order 
to construct generalized scheme based on the evolution of polynomials of higher order. In this 
section we explore the case of cubic polynomial 

p{X) = X^- 3Pi X'^ + 2P2X -P^ = 0, (5.1) 

relations between coefficients and eigenvalues are given by 

3Pi = 5i g2 + gs, 2P2 = qiq2 + q2q3 + qsqi, P^ = gig2?3- (5.2) 

We shall restrict ourselves only with the case when coefficients are represented by real numbers 
and let us assume that p{X) is irreducible. By replacing X with X — Y + Pi we come to 
invariant polynomial 

Y^ + R2Y -Ro^O, (5.3) 
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where 

(a) R2 = 2P2 - 3Pl^ ih) Pf + i?2 Pi + Po = P^- (5.4) 

The eigenvalues and the coefficients of this polynomial are invariants with respect to simulta- 
neous translations of the eigenvalues of Eq. (5.1). Obviously this statement is a consequence 
of the formula Y = X — Pi from which it follows 

3yi = 62 - 63, 3y2 = 63 - ei, 3y = ei - 62, (5.5) 

where 

ei = ?3 - ^2, 62 = ?i - ?3, 63 = ?2 - gi. 

The evolution equations remaining constants coefficients Pq, R2 are obtained directly from 
Eqs.(5.4). Differentiating these equations we get 

dPi {3P^ + P2) = 2P2 dPi = dP^, (5.6) 

dP2 = SdPi. 

Evolution equations for the eigenvalues, obviously, have to has the following form 

-^Qk^A, k= 1,2,3, 
as 

where A is some function same for all eigenvalues. In order to construct some dynamics in the 
quahty of A we must take A — P. Then 

The cubic polynomial is an object of special interest, because this polynomial is closely 
related with the classical elliptic functions. The solutions of the evolution equations for the 
eigenvalues and the coefficients of the cubic polynomial can be represented via elliptic Jacobi 
and Weierstrass functions, correspondingly [1]. On making use (5.7) from (5.4b) we come to 
the following differential equation 

dP 

Pl + R,P, + R, = {-lf. (5.8) 

ds 



Write this equation in the following designations 

Ap{zf - g2p{z) - 5/3, (5.9) 



dz 

where z — 4s, g2 — — 4Pi, — — 4Po and p(2s) = -Pi(s). The integral formula for p{z) is 
given by 

z^ {ix^ - g2X - g^y^/^dx. (5.10) 

The functions p{z) and p'{z) are Weierstrass elliptic functions with periods 2a;i, 2a;2. Define 
cus — —cui — CU2- Then the values 

p{ui), p{u2), p{uJ3) 

are the roots of cubic equation 

ix^-g2X-gs^0. (5.11) 
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Introduce new variables pk, k = 1,2,3 where pi — Qk- Then P — piP2P3- For p^ evolution 
equations are derived from (5.7): 

dpi dp2 dps 

-T = ^'SPS, -r = PlPS: -r ^ (5.12) 

as as as 

Solutions of these equations presented by quotients of Jacobi elliptic functions. Let sn{u), cn{u), dn{u) 
be a set of Jacobi elliptic functions. Define the following quotients of these functions (in Glaisher 
notations [1]) 

1 cn dn 

ns — , cs — , as — , 

sn sn sn 

which obey the following differential equations 

—ns(u) = cs(u)ds(u), —cs(u) = ns(u)ds(u), —ds(u) = cs(u)ns(u), (5.13) 
du du du 

with 

ns^ — cs^ = 1, ds^ — cs'^ = 1 — k. 
Prom Eqs.(4.12) it follows that the values 

ei = - ?2, 62 = gi - Qs: 

are constants of motion. Hence solutions of Eqs.(4.12) are presented via Jacobi elliptic functions 
as follows 

qs = eins^{u,k), q2 = eics'^{u,k), qi = eids'^{u,k), k=l . (5-14) 

ei 

The final result let us represent by the following 
Statement: 

// squared roots of the eigenvalues of the cubic equation obey equations for Jacobi elliptic 
functions, then the evolution of the coefficients are governed by equations for Weierstrass elliptic 
functions. 

For cubic equation there exist, also, another possibility to express its eigenvalues, namely, 
via trigonometric functions. To the formulae given below we come virtue of formulae (3.25) 
derived in Ref.[8]. 

The eigenvalues of the reduced polynomial are given by formulae: 

yi — — ^^3i?i cos^, y2 — ^\fsRi{cos9 + V3sm9), y^ — ^'\fsRi{cos9 — VSsinO). (5.15) 
It is easy to verify that 

yi + y2 + y3 = 0, yiy2 + y2y3 + y3yi = -Ri- 

Equation for the last coefficient, Rq, leads to trigonometric equation for cos(^): 

-Ro ^ ymys 

97 R 

= COS e(cos^ ^ - 3 sin^ 9)^4 cos^ ^ - 3 cos ^ ° 



2 ivmr 

By using the trigonometric formula 

cos (3^) = 4cos^e - 3cose 
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we come to simple trigonometric equation: 



cospO) = 5^^. (6.16) 



Since we restrict ourselves only with real solutions the following inequality has to be true 

(^r < i^r- (5.17) 



The algorithm elaborated in the previous section in the case of cubic equation (5.1) is 
simplified as follows. With respect to x — as an evolution parameter evolution equations 
(5.6)-(5.7) have a form 

2P2dPi = 2P2dx, 2P2dP2 = SPidx, (5.18) 

with the initial data Pi{x = P^) = Pi, P2{x = P^) = P2. At the final stage when a; = P^ = 

one finds Pi(P — 0) — Pi, P2(P — 0) — P2 satisfying relationships 

(1) (1) (1) (1) 
R2 = 2P2 - 3P2, P3 + R^P^ + R^ = 0. 

Thus the new equation possesses with same invariants as the original one, 

(1) o 

^3(0) -3PiX2(0) + 2P2X(0) =0. (5.19) 

Equations (5.1) and (5.19) possess with congruent eigenvalues, however Eq.(5.19) has one trivial 
root. Therefore the problem is reduced to the solution of quadratic equation 

X\0) - 3PiX^(0) + 2P2X(0) = 0. (5.20) 

Prom solutions of eq.(5.20) to the solutions of Eq.(5.1) one comes simply by the set of transla- 
tions: 

(1) ^ ^ (1) ^ ^ (1) 

gi = Pi-Pi, g2 = 52(0) + Pi -Pi, g3 = ?3(0) + Pi-Pi. 

The dynamics related with the cubic polynomial, evidently, is characterized with two first 
constants of motion P21 -^o which do not depend of potential field. In quality of evolution 
equations we have to use Eqs.(5.6)-(5.7) written with respect to parameter of evolution defined 
in (4.6): 

^2 = ^P, = 3dPi, -^Pi = P (5.21) 

as as as 

Evolution equations for the eigenvalues, correspondingly, are given by 

^pI = P. A; = 1,2, 3. (5.22) 
Prom these equations it follows that the constants of motion are given by formulae 

Ml = p2 _ pI M2=pI- pI, Ms=pI- pI (5.23) 
In order to include into this scheme the potential field V{r) we have to use equation (4.6): 

— = (U-n), \] = VVir), Pn = P. (5.24) 
or 
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Furthermore, the set of evolution equations has to be completed with the interrelation between 
momentum and velocity with respect to time-like parameter: 

P ^ ^. (5.25) 

In these designations evolution equations (5.21) are transformed into the following dynamic 
equations 

^ _VV(r) P2, = -3(P • V)V(r)Pi, ^Pi = -(P • V)V(r). (5.26) 

In the case of stationary potential, beside the first constant R2, Ro, the equations imply another 
constant of motion, the energy 

£^ = + V. (5.27) 
On making use of expression for Pi from (5.27) in the second equation of Eqs.(6.6), we get 

£2 = P2 + 3£iV - ^V^ = 1{R2 + 3^1'). (5.28) 

This is second expression for the energy. Next, express from this formula P2 and substitute 
into first equation of (5.26). This leads to Newtonian equation 

A.p^-^W{r,S^), (5.29) 

with effective potential 

2W{r,Si) = S2V{r) - SSiV^r) + V^{r). (5.30) 
From Eq.(5.29) one may find formula for the total energy 

S^^P^ + W^^{R^ + + Sf). (5.31) 

6 Generalized dynamics related with evolution of n- 
degree polynomial 

Now we have accumulated enough experience in order to be able to build the general form 
of the dynamics related with n-degree polynomials. Let us start with evolution equations 
(2.8) written for n-dcgree polynomial (2.1). Dynamic equations with respect to some time-like 
evolution parameter s describing a motion inside stationary potential field V{r) are formulated 
in the following form 

^ = U P„_i, (6.1a) 

dPk 



(U-P)Pfc_i(n- A; + 2), A; = 2, . . . , n - 1, (6.16) 
= (U-P). (6.1c) 



dr 

dPi 



dr 

The set of evolution equations has to be completed with the interrelation between momentum 
and velocity given by Eq.(5.25). This system of dynamic equations remain invariant the coef- 
ficients RQ,Rk,k — l,...,n — 1. The motion of a physical system obeying to these equations 
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possesses with the set of inner and outer momenta. Evolution of the set of outer-momenta 
{P^, Pk, A; = 1, .... n — 1 } are given by Eqs.(6.1), whereas the evolution of the inner-momenta 
{pk, k = 1, ...,n} is described by 

^ = (U-P), k^l,...,n. (6.2) 

The first integrals of this system are given by 

Mik^pl-pl i^k. (6.3). 

Prom Eq.(6.2c) we find the first constant of integration, the energy Si = Pi + V, or Pi = 
£i — V. By substituting Pi into the next, k — 2-th equation of (6.2b) we find the other constant 
of integration (the second expression for energy): 

S2 = P2 + nEi y -i^V". 

By continuing this process, namely, by substituting P2 from the last expression into the next 
equation with A; = 3 we shall find the third constant of integration: 

£3 = _ („ _ V + ^^E, V- - '^^V\ 

Continue this process up till {n — l)-th stage. Then, in the {n — l)-th stage we shall obtain 
the expression for Pn-i- By introducing this expression into Eq.(6.2a), we come to Newtonian 
equation 

Ap = _vw^(r,£i), (6.4) 
ar 

where the effective potential is given by the following series 

1 " 1 
W{r,£i) = -Y,{k + l)£^_kV\r){-lf^\ with £0 = (6.5) 

Prom Eq.(6.4) one may find formula for the total energy 

1 1 ""-^ 

£^-P^ + W^ -{Ro + E RkSf + ^r)- (6.6) 

^ ^ k=2 

Concluding remarks 

This is a prodigious fact that the evolutions equations elaborated for polynomials serve as 
dynamic equations for the generalized dynamics of high-energy particles. Notice, the Algorithm 
of calculation of eigenvalues of the polynomials elaborated in this paper is distinct of the 
numerical methods of calculations which principally arc based on iteration process with initial 
sampling, or tentative, data of roots, so that effectiveness of these algorithms depends of the 
initial data. The iteration process based on the present Algorithm uses in the quality of initial 
data the coefficients of the original polynomial. The method elaborated here can be considered 
also as a theory of functional connection between coefficients and eigenvalues of the polynomial 
expressed as one-valued Weierstrass and Jacobi hyper-elliptic functions. Evidently, the present 
method without any principal difficulties can be continued to the case of polynomials defined 
over the field of complex numbers. 

We have restrict our attention only on dynamic equations given in one-dimensional coor- 
dinate space. On the theory of the generalized dynamics in AD space with physical units one 
may consult in Refs.[15] and [16] and the references therein. 
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